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0> ' Abstract: By solving the Einstein equations of the graviton couphng with a real 

OO ' 

^ ■ scalar dilaton field, we establish a general framework to self-consistently solve the 

geometric background with black-hole for any given phenomenological holographic 
models. In this framwork, we solve the black-hole background, the corresponding 
dilaon field and the dilaton potential for the deformed AdSs model with a posi- 
tive/negative quadratic correction. We systematically investigate the thermodynam- 
ical properties of the deformed AdSs model with a positive and negative quadratic 
X : correction, respectively, and compare with lattice QCD on the results of the equation 

. of state, the heavy quark potential, the Polyakov loop and the spatial Wilson loop. 

We find that the bulk thermodynamical properties are not sensitive to the sign of 
the quadratic correction, and the results of both deformed holographic QCD models 
agree well with lattice QCD result for pure SU(3) gauge theory. However, the results 
from loop operators favor a positive quadratic correction, which agree well with lat- 
tice QCD result. Especially, the result from the Polyakov loop excludes the model 
with a negative quadratic correction in the warp factor of AdSs. 
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1. Introduction 



In recent decade, the discovery of the anti-de Sitter/conformal field theory (AdS/CFT) 
correspondence and the conjecture of the gravity/gauge duahty has been widely 
used to understand strongly coupled dense and hot quark matter [0, |, |, |], ^ ^ |], 

1,0,0,0,0 0. 

One of the most successful examples of using thermal A/" = 4 super- Yang-Mills 
theory (SYM) in quantum chromo dynamics (QCD) at finite temperature is the small 
value of shear viscosity over entropy density rj/s = l/An [Q, which is very close to 
the value used to fit the RHIC data of elliptic fiow V2 10, It is now believed 
that the system created at RHIC is a strongly coupled quark-gluon plasma (sQGP) 
and behaves like a nearly "perfect" fiuid |0, IE . However, lattice QCD results show 
that the bulk viscosity over entropy density ratio (/s rises dramatically up to the 
order of 1.0 near the critical temperature [jl9|, ^ The sharp peak of bulk 



viscosity at Tc has also been observed in QCD effective models ||2^, The large 

bulk viscosity near phase transition is related to the non-conformal equation of state 
[ p5| , p6| , e.g, the square of sound velocity is around 0.07, which is much smaller than 



the conformal value 1/3. In the SYM plasma, the square of sound velocity is always 
1/3 and the bulk viscosity ( always vanishes at all temperatures. 

These facts demonstrate that strongly coupled quark gluon plasma near Tc is 
not conformal invariant. In order to mimic the QCD equation of state, much effort 
has been put to find the gravity dual of gauge theories which break the conformal 
symmetry, e.g, p, |], ^ ||, where a real scalar dilaton field background has been 
introduced to couple with the graviton. Refs. 0] and p, ^ have used different dilaton 
potentials as input, while Ref. has used QCD /3-function as input. 

The AdS/CFT approach has also been widely applied to describe non-perturbative 
phenomenology in the vacuum, where conformal symmetry is also needed to be bro- 
ken 



m 40, 41, 42, 43, 44]. The most 



>8 29, 30, 31, 32 



33| g4|, |35|, ^ |37|, |S 

economic way of breaking the conformal symmetry is to add a proper deformed warp 
factor in front of the AdSs metric structure, which can capture the main features of 
non-perturbative QCD physics. For example, Andreev-Zakharov proposed a positive 

2 

quadratic correction, e'^^ with z the fifth dimension coordinate and c > 0, in the 
deformed warp factor of AdSs geometry, which can help to realize the linear heavy 
quark potential The linear heavy quark potential can also be obtained by in- 
troducing other deformed warp factors, e.g. the deformed warp factor which mimics 
the QCD running coupling and the logarithmic correction with an explicit IR 
cutoff log ^'^^^ p4| . To produce the linear Regge behavior of the hadron excitations. 



Karch-Katz-Son-Stephanov |34| proposed the soft-wall AdSs model or KKSS model 
by introducing a quadratic dilaton background in the 5D meson action, whose effect 

2 

in some sense looks like introducing a negative quadratic correction, e""^^ , in the 
warp factor of the AdSs geometry. However, it is worthy of mentioning that the 
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model with a quadratic correction in the warp factor of the metric is not equiva- 
lent to the model with a quadratic correction in the dilaton background. A positive 
quadratic correction e'^^ in the dilaton background of the 5D hadron action has also 
been used to investigate hadron spectra 0, however, higher spin excitations in 
this background will get imaginary mass ^ . 

From the above studies of constructing the holographic QCD (hQCD) models for 
describing the heavy quark potential and the light hadron spectra, we have observed 
that a quadratic background correction is related to the confinement property, i.e. 
the linear quark anti-quark potential and the linear Regge behavior. Furthermore, we 
have observed that to produce the linear heavy quark potential, the positive quadratic 
correction has been introduced in the deformed warp factor, while to produce the 
linear Regge behavior, the negative quadratic correction is introduced in the dilaton 
background in the 5D hadron action. 

It is natural to ask the connection between the quadratic correction in the met- 
ric and dilaton background, and which sign is consistent with real QCD data. If 
we assume that there is no direct connection between the dilaton field and metric 
background, these two types of models cannot be contrasted. However, the metric 
structure of the hQCD model and its corresponding dilaton background can be self- 
consistently solved from the Einstein equation in the non-critical string framework 



or the 5D Einstein dilaton system, e.g, as described in Ref. In this frame- 

work, the connection between the metric structure and the dilaton background can 
be established. 

In Ref.[^, we have numerically solved the dilaton background at zero temper- 
ature for the deformed AdS5 model with positive and negative quadratic correction, 
respectively. In this paper, we extend our previous study in the non-critical string 
framework at zero temperature to finite temperature, and self-consistently solve the 
black-hole background, dilaton field and the dilaton potential. In order to further 
pin down the sign of the quadratic correction in the hQCD model, we systematically 
investigate the thermodynamic properties of these two deformed AdSs models, in- 
cluding the equation of state, the heavy quark potential, the spatial string tension, 
and the Polyakov loop at finite temperature. By comparing with lattice results, we 
find that the equation of state is not sensitive to the sign of the quadratic correction. 
However, the results from loop operators favor a positive quadratic correction. Espe- 
cially, the result from the Polyakov loop excludes the model with negative quadratic 
correction in the warp factor of AdSs. 

The paper is organized as follows. In Sec.|], we derive the general formulae for 
equations of motion in classical 5D gravity-dilaton system, and solve the black-hole 
background, the dilaton field and its corresponding dilaton potential for any given 
hQCD model. In SecH, we self-consistently solve the dual black-hole background, the 
dilaton field and its dilaton potential for the hQCD models with positive and negative 
quadratic correction in the deformed warp factor of AdSs metric background. In 
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Sec.^, we study the phase transition and the equation of state, including the entropy 
density, the pressure density, the energy density and the square of sound velocity for 
the two hQCD models with quadratic corrections. In Sec.^, we calculate the heavy 
quark potential, the Polyakov loop and the spatial Wilson loop at finite temperature 
and compare our results with lattice data. The summary and discussion is given 
in Sec.^. We also list several special exact solutions of the general graviton-dilaton 
system in Appendix ^ and the method of extracting the spatial string tension in 
Appendix 0. 



2. 5D Einstein dilaton system with black-hole background 



2.1 Minimal non-critical string framework 

For low energy IIA/IIB supergravity, the classical action in the string frame has the 
following form 



SiiA/IIB — Sns + + Scs + Sf^ 



ermion- 



(2.1) 



Where Sf ermion is the action from the fermionic part, and S^s^ Sr, Scs are bosonic 
part actions for Neveu-Schwarz (NS) sector of graviton-dilton coupling, Ramond (R) 
kinetic term for the form field, Chern-Simons term, respectively. The NSNS sector 
action takes the form of 



S 



NS 



2 ' ^' 



(2.2) 



The RR part action and Chern-Simons term take different forms for Type-IIA and 
JIB supergravity, respectively: 



Scs\iiA 
Sr\iib 
Scs\iiB 



1 



9^2 



9P 



^2 A F4 A F4, 



(2.3) 



Where kiQ the 10-dimension Newton coupling constant, Gs the metric in the string 
frame, R, (p are scalar curvature and the dilaton field, respectively. B is anti- 
symmetric 2-form tensor field, C is R-R form field, H, F are field strengths of B,C, 
and F = e~'^F. ^2,^4 only appear in IIA and ^1,^3,^5 only appear in IIB super- 
gravity. 
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We are interested in building a 5-dimensional effective gravity dual to a 4- 
dimensional non-Abelian gauge theory. We follow the argument in Ref . to derive 
the minimal 5D effective graviton-dilaton system from lOD supergravity Eq.( |2.1| ). 

1) For pure Yang-Mills gauge theory, there is no fermionic degrees of freedom, 
therefore, we neglect the Sjermion in Eq.(|2.1|) for simplicity. 

2) In the NSNS sector, one should note that we ignore the effect of = dB2, 
which is normally taken to be zero because it induces noncommutative effect in the 
string theory. This leads to the vanishing of the Chern-Simons action, i.e. Scs = 0. 

3) The RR part includes six fields strengths ^0,^1,^2,^3,^4,^5, where only 
Fq, Fi, F2 are independent due to the Poincare dual relationship. Fq couples to bulk 
instantons, and its dual F^ Fq generates a four-form field stems from branes 
related to the U{Nc) gauge group. Fi = dCo originates from axion field Cq = a 
which is dual to Tr[F A F] in the gauge theory. F4 ~ Fi and its associated three- 
form C3 couples to domain walls. F2 generates a vector and couples minimally to the 
baryon density, and it breaks Lorentz invariance in the QCD vacuum. Therefore, F2 
and its dual F3 can be safely neglected. 

From the above analysis, the relevant low-lying fields, which are expected to be 
dual to the lowest dimension operators in the gauge theory, are the graviton g^^, the 
dilaton 0, and the RR axion a. One can put the solution of F5 into Eq.( p.l|) and 
then integrate out the inner five compact space-time M5. Finally, one can obtain the 
minimal non-critical 5D effective gravity action as following: 



S5D = j (f'x^/^e-^'f' {R + ^d^4>d^<\> - Vsm 



(2.4) 



Where is 5D Newton constant, and ^(0) are the 5D metric and dilaton 
potential in the string frame, respectively. 

In the minimal non-critical string framework, the metric structure, the dilaton 
field and the dilaton potential should be solved self-consistently from the Einstein 
equations. One can use the dilton potential as input to solve the metric structure and 
the dilaton field. There are several different choices for the dilaton potential from 
several groups. The scalar field or dilaton field encodes the running of the Yang- 
Mills gauge theory's coupling A. In the framework of Giirsoy -Kiritsis-Nitti (GKN) 
[pl , ^ the renormalized dilaton field has been defined as A = e"^, and the dilaton 
potential can be solved from the input of QCD running coupling function, i.e, /3(A)- 
function. In Ref. [§, in order to mimic the QCD equation of state, the authors 
chooses the dilaton potential as ^(0) = where L is the radius of the 

AdSg. 

However, in the graviton-dilaton coupling system, for the given dilaton potential 
in H, ^ the solution of the metric structure is very complicated. Moreover, 
for finite temperature case, the temperature dependence of the dilaton potential is 
not obvious. In Ref. [HI , instead of input the dilaton potential, we solve the diltaon 
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potential for given metric structure of deformed AdSs model, which can capture some 
features of QCD phenomenology at zero temperature. In this work, we will extend 
the minimal non-critical string framework to finite temperature, and investigate the 
phase transitions, the equation of state and the temperature dependence of loop 
operators. 

2.2 Graviton-dilaton system with dual black-hole for hQCD model 

The finite temperature dynamics of gauge theories, has a natural holographic coun- 
terpart in the thermodynamics of black-holes on the gravity side. Our holographic 
model is defined in the string frame, it is convenient to calculate the vacuum expec- 
tation value of the loop operator in the string frame, however it is more convenient 
to work out the gravity solution and study equation of state in the Einstein frame. 
Therefore, for later use, we will firstly derive the relation between the string frame 
and the Einstein frame. 

If the metric in the Einstein frame g^^, and its corresponding metric structure in 
the string frame g^^, are connected by the scaling transformation 

4. = (2-5) 

then the scalar curvature in the Einstein frame and string frame has the following 
exact relationship: 

g-2a^5 ^ _ (£) _ _ 2)du,fld^n + 2{D - (2.6) 



with D the dimension, and is defined by / dn^/—g^d'^, which is useful to 

derive the exact relation between two actions in string frame or Einstein frame. 
In the case oi D — 5, we have 

J y^e'^^i?^ = J y3^e('"-2)^ [i?^ - Ud^Wn + 8(m - 3)V^n] , (2.7) 

J V^e"^^ {g's'^d.ndM) = J v^e(-^)^ {e'^'g^/d.Qd^Q) . (2.8) 
The dilaton potentials in two different frames own the following relationship 

Vs{n) = VE{n)e^''. (2.9) 
By setting m = 3 and fl = — 10, we have 



Vs^VEe—, (2.10) 



and the following relation: 



/ 
/ 



(2.11) 
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Therefore, Eg . (12 .41 ) in the Einstein frame becomes: 

S,D = / d'x^/^ (^R - ^9^09^0 - VEi(l>)^ (2.12) 

Adding the black-hole background to the holographic QCD model constructed 
from vacuum properties, in the string frame we have 

dsl = -Ar- -f{z)dt^ + TTT + dx'dx' , (2.13) 



with L the radius of AdSs. The metric in the string frame is useful to calculate the 
loop operator. To derive the Einstein equations and to study the thermodynamical 
properties of equation of state, we transform it to the Einstein frame 



dsl = \^-f{z)df + j^ + dx'dx'j . (2.14) 

The general Einstein equations from the action ( |2.12| ) takes the form of: 

E,u + (^9^</'5'^0 + VEi<P)^ - \d^4>QA = 0. (2.15) 



By using Eqs.( p.l4D and ( p.l2| ), we can derive the following nontrivial Einstein equa- 



tions in {t,t), {z,z) and {xi,Xi) spaces, respectively: 

A'liz) + Kiz) ( §^ - U'iz) - 2 + A'^z)) - f{z) ^ 



2f{z) 3^ ' ' z ''7 ' ' ' K^fiz) 2zf{z) 



_|ni + l + £;M*(£))e-...-^.0 (2.17) 

Arf{z) z^ 4:Z^f{z) 



+/(.)(6A;(.)-^-^-MPfM + ^'+6^'.(.f 

- 40"W + 4,#.'(j)2 + M£) ) = (2.18) 

Z 

One should notice that the above three equations are not independent. There are 
only two independent functions, one left equation is to check the consistence of the 
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solutions. To simplify the above equations, one can use the following two equations 
without dilaton potential Ve: 

4>"{z) - {2K{z) - - ^ + - ^ = 0, 

f"{z) + (^A^z) - 20'(^) - f'{z) = 0. (2.19) 

The EOM of the dilaton field is given as following 

^3 ^^ ^ d.<t>j - O^Ve = 0. (2.20) 

For any holographic QCD model with given metric structure As{z) in the string 
frame, we can derive the general solutions to the Einstein equations, which take the 
following form: 

- e^M-) 3AJz) 
^0 + 01 / —dx + 







x2 2 
3 r e"^^^-^ j^y^e-^^'^y^A^yfdy 



2 Jo 

f{z) = fo + fJ [ a;3e2'^(^)-3^=(^) dx 



dx. 



Ve(0) = — 



\J0 



L2 

z'fiz) - 4/(z) {3z'A':{z) - 2z'<p"{z) + z^<P\zf + 3)) . (2.21) 

Where 0o;0i;/o)/i are constants of integration. In Appendix 0, we list several 
simple exact solutions of the Einstein equations by using Eq. ( ^.211 ). 

3. Dual black-hole solution for the hQCD model with quadratic 
correction 

From the experiences of constructing holographic QCD models for describing the 
heavy quark potential and the light hadron spectra, we have learnt that a quadratic 
background correction is related to the confinement property, i.e. the linear quark 
anti-quark potential and the linear Regge behavior. A positive quadratic correction, 
e'^^^ with c > 0, in the deformed warp factor of AdSs can help to realize the linear 
heavy quark potential |^|. A quadratic dilaton background in the 5D meson action. 



whose effect in some sense looks like introducing a negative quadratic correction, 
e~'^^ , in the warp factor of the AdSs geometry, is helpful to realize the linear Regge 
behavior of hadron excitations 
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Therefore, we introduce the following holographic QCD model with positive and 
negative quadratic correction in the deformed warp factor of AdSs in Eg. (|2.13|) , i.e. 

we take 

A,{z) = ck'^z^ (3.1) 

with c = ± indicating the positive and negative quadratic correction, respectively. 
The main purpose in this section is to solve the dual black-hole background, corre- 
sponding dilaton field, and dilaton potential in the Graviton-dilaton system for the 
given hQCD model Eq.(|3.1|). To get the solution of equations Eq.( |2.16| )-Eq.( p.20| ), 



we impose the asymptotic AdS^ condition /(O) = 1 near the UV boundary z ~ 0, 
and require 0, / to be finite z = 0, Zh with Zh the black-hole horizon. Fortunately, 
we find that the solution of the black-hole background takes the form of, 

f{z) = l-f^ j'j exp 0ca;2 {H,{x/k) - 1)^ dx, (3.2) 

with 

fc = — ^ ' (3-3) 

/q x^exp (fcx^ {Hc{x/k) - 1)) dx 

and 

Hc{z)= 2F2(^l,l;2,|2cA;V^. (3.4) 

It is worthy of mentioning that in Ref . |TD|, |TT], , the authors also investigated 
the thermodynamic properties of the deformed AdSs model with positive quadratic 
correction in the warp factor, however, there the black-hole background was not 
solved self-consistently from the graviton-dilaton system, but just taken from the 
AdSs Schwarz blak-hole (AdS-SW BH) background 

fiz) = 1 - zV4, (3.5) 

with zt = 1/ttT. 

The solutions for the dilaton field and the dilaton potential take the forms as 
following: 

(j){z) = lck^z\l + H,{z)), (3.6) 



64cA;V - 38Ak^z^ + 12V2^e^'''' kz{ -7 + 20ck^z^)Erf,{V2kz) 
-27716^'''''' Erf, (^V2kz^ - 

ofhlck2z^~l+H4z))r.3 3 

^ ^ Ukz - 16cPz3 + 3V2^e^''' ' ErfJV2kz)] 

(3.7) 
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We have used to represent the solution of the dilaton potential for c = ±. 
Where the function 2-^2 ({oi, .., flp}; {&i, feq} ; -z) is a generalized hypergeometric 
function which has series of J^'kLo ("^Ofc • • • i'^p)k I (^i)fc ■ ■ • (^9)^ I ^' /o is related 
to the integral constant which can be expressed by the position of the horizon of the 
black hole solution. The function Erfc[z] appeared is also a special function class 
named by error function which is defined as integral form Erf^[z] = Erf[z] = 
^ /; e-'"dt and Erf_[z] = Erfi[z] = 

Conformal invariance in the UV can be obtained when ~ at the UV boundary 
2 — )■ 0. One can expand (f){z) at UV boundary z ~ 0, 

0(^-^0) _ + __ + .... (3.8) 



The behavior shown in Eq. ( |3.8| ) is consistent with the requirement of the asymptotic 
AdSs near the ultraviolet boundary. 

Through the AdS/CFT dictionary, for any dilaton field $, we have 

lim V(<|.) = + ^A(A - 4)$^ + 0($^). (3.9) 
By using the below relationship, 

^2,,riT.^ _ dr d f dr dV{z)\ _ ^^2^2 



one can easily get the conformal dimension of the dilaton field near the UV boundary. 
In Eq. ( p.9|) , A is defined as A(A — 4) = M|L^, which is constrained by the 
Breitenlohner-Freedman (BP) bound 2 < A < 4. In our case, $ = 

Ml = --i, (3.11) 

for c = ±. Therefore, the conformal dimension of the dilaton field is A = 2 in our 
holographic models with c = ±, which satisfies the BP bound but does not correspond 
to any local, gauge invariant operator in 4D QCD. Although there have been some 
discussions in recent years of the possible relevance of a dimension two condensate 



in the form of a gluon mass term [J7 , it is not clear whether we can associate with 
dimension-2 gluon condensate, because the AdS/CFT correspondence requires that 
bulk fields should be dual to gauge-invariant local operators. 



4. Equation of state for hQCD model with quadratic correc- 
tion 

A black-hole solution with a regular horizon is characterized by the existence of a 
surface z = Zh, where f{zh) = 0. The Euclidean version of the solution is defined 
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only for < 2; < z/i, in order to avoid the conical singularity, the periodicity of the 
Euclidean time can be fixed by 



r r + 



\f'i^H)\- 

This determines the temperature of the solution as 

\f'{^k)\ 



(4.1) 



T 



An 



(4.2) 



From Eq. ( |3.2| ), one can easily read out the relation between the temperature 
and position of the black hole horizon. 



^_ k'zlexp{lick'zlH^iz,)-cezl)) 
4.1 The would-be critical temperature 



(4.3) 



According to Eq.( [4.3D , we cannot get the analytical expression to describe the relation 
between the temperature and horizon. The numerical results for the temperature as 
a function of the horizon Zh for c = + and c = — are shown in Fig.|l] (a) and (b), 
respectively. Here we have taken k = 0.43GeV. The dashed lines are results from 
the pure AdSs Schwarz black- hole in Eq.( p.5|) . 

T(GeV) 



0.6 



0.4 




AdS-SWBH 



;i(fm) 




AdS-SWBH 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

(a) 



;ilfm) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

(b) 



Figure 1: The temperature T as a function of the black-hole horizon with k = OASGeV 
for c = -|- in (a) and c = — in (b), respectively. The dashed Hues are for pure AdSs Schwarz 
black-hole. 

From Fig. |l], it is noticed that for pure AdSs Schwarz black-hole, the temperature 
monotonically decreases with the increase of the horizon [|Ty, |rT|, ^ . If one solves the 



-11- 



dual black-hole background self-consistently, one can find that for both hQCD models 
with positive and negative quadratic correction, there is a minimal temperature 
Tjnin at certain black- hole horizon This is similar to the case for the confining 
theory (at zero temperature) discussed in Ref. 0. For T < Tmm, there are no 
black-hole solutions. For T > Tmm, there are two branches of black- hole solutions. 
When Zh < z^, the temperature increases with the decreasing of Zh, which means 
that the temperature increases when the horizon moves close to UV, this phase 
is thermodynamically stable. When Zh > z^, the temperature increases with the 
increase of Zh, which means that the temperature becomes higher and higher when 
the horizon moves to IR. This indicates that the solution for the branch Zh > z^ is 
unstable and thus not physical. 

In order to determine the critical temperature, we have to compare the free 
energy difference between the stable black hole solution and the thermal gas. The 
thermal gas is solved by setting /o(-2) = 1 0, From the AdS/CFT conjecture, 
the gravity side is weakly coupled, so we can use the semi-classical approximation 
and just deal with the on-shell action. Following Ref. [53] and Appendix C in Ref. 



the regularized total on-shell action in Euclidean space can be decomposed as 
three parts: 

tR tR I tR I tR _L _L jcount ( A A\ 

— ^ ^GH ^count — "I" "I" , {^■^) 

with the on-shell Einstein action, Iq^j the Gibbons-Hawking term and I^^nt 
counter term from the renormalization scheme. P and are the contributions at 
UV cut-off z = e and IR cut-off z = zir, respectively. Correspondingly, the free 
energy density is defined as 

jr ^ ^jR = jre^ jzIR + jzcount^ ^4 5^) 

with V3 the volume of the three-space. 

To compare the free-energy density of the black-hole case and the thermal gas 



case, one has to impose the following conditions ||49|| : 



0o(z) 1= <P{z) |„ (4.6) 

M^) 1= f5h{z)^J{z)l, (4.7) 

V^hl{z) 1= V^h\z) U . (4.8) 

Where (/)o, &05 /3, V3 are for thermal gas with 60 = ^e^^°*^^) and Aeq = Ago — |</)o, and 
(j), b, /5, V3 are for the black-hole with b{z) = ^e^^^^\ 

The free energy density for the black-hole has the form of 

J'bh = nn + ^bT' (4.9) 

nu = (3b{erf{e)b'{e) + ^^(e) ^(e)) , (4.10) 
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where we have defined = l/(167rG5). It is noticed that for black-hole case, zjr 
is normally set at the horizon Zh, therefore vanishes due to f{zh) = 0. At the 
limit e — >■ 0, we have 

^hn = -^i^^ - l^h'L'M' - h^L^M^f^'^ + o{e% (4.11) 

For detailed derivation of Eq.([4.11[), please refer to Appendix O. 



The free energy density for the thermal gas takes the form of 



^TG — ^TG + -^TG + ^TG^^i (4-12) 

T'rl = 2M^U{ziR)hl{zm%{zm). (4.14) 



e rv/o e V 2 y 



At the limit e — )■ 0, from the results derived in Appendix the free energy density 
of the thermal gas at UV cut-off takes the form of 

r^G = - (^Ik'LV - Ik'LVf^^ + o{e% (4.15) 

Where we have used /o(^) = 1 in the derivation. 

The divergent behaviors of the thermal gas and black-hole are the same, i.e, the 
counter terms J-'^^"'^ and J-j.^^"* would have the same contribution in the two phases 
thus J-"^"^"* — J^^°q"'^ would disappear in the free energy difference. Therefore, the 
final expression of the free energy difference has the form of 

AT = Tbh-J'tg= -\k^M^L^f^^ -2MX{zm)h',{zjR). (4.16) 

It is worthy of mentioning that > 0, therefore the contribution from the first 
term is always negative. Furthermore, from Eq.( B.15 ), we can observe that the first 



term contribution is exactly —\Ts with s the entropy density, and the free energy 
difference can be written as: 

AJ- = -^Ts - 2M%l{ziR)h',{zjR). (4.17) 

Because h^^z) is monotonically decreasing with z so &o(-2) < 0, the second term con- 
tribution to the free energy difference is always positive. The expression of Eq.( |4.17l ) 
is similar to Eq.(3.25) in Ref. [Q. We can observe that the second term contribution 
m Eq.(|]T^ plays the same role of the gluon-condensate in Ref.fg]. 

The second term contribution in Eq. ([4.17| ) is zjr dependent. If we follow Ref.0 



and choose zir at the good singularity point, i.e, zir — )■ oo, where bo{ziR) — 0, then 
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the second term contribution in Eq.( 4.16 ), i.e, the zjr related term vanishes and we 
have 



-Ts < 0. 



(4.18) 



This indicates that the black-hole phase is more stable than the thermal gas. The 
free energy difference for the case of zjr — )■ oo as a function of T/Tmin is shown in Fig. 
1^. The solid lines indicate the stable black-hole solution, and the dashed lines are for 
the unstable black-hole phase. It is seen that AJ^ jumps to a negative value at T^in, 
then decreases with the increase of the temperature. In this case, we observe a zeroth 
order phase transition at T^m due to the free energy discontinuity. By choosing the 
parameters k = 0.43GeV and G^/L^ = 1.26, the value of — AJ-'/T^ approaches 1.5 
at high temperature, which is similar to the lattice data of the pressure density for 
pure SU(3) gauge theory . 




(a) 




Figure 2: The free energy difference of Eq.( [4.18D scaled by as a function of T/Tmin 
with k = 0A3GeV and G^/L^ = 1.26 for c = + in (a) and c = — in (b), respectively. The 
dashed lines are for unstable black-hole solution. 



In principle, the critical temperature from Eg. ( [4. 17] ) should be determined 
by finite value of zjr, which sets the basic scale of the theory and should be fixed 
by experimental or lattice data. For example, for zjr = 2.34GeV~^, the free-energy 
difference will cross the zero-energy axis at Tmm so that Tc = Tmin and the phase 
transition will be of first order. If we choose other smaller value of zir, e.g, zir = 
2GeV-\ the free energy density difference will cross the zero-energy axis at > T^m 
and the phase transition will be also of 1st order. 



For our numerical calculations, we will firstly choose T, 



as the "would-be" 

critical temperature, which means that we choose a finite value of zjr = 2.34GeV~^ 



min 
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and make 5 J-" = at Tmm- The exact thermal gas solution in the region T < T^m is 
unknown for us, therefore, we only focus on the temperature region T > Tmin for the 
deconfined quark gluon plasma phase. In the following we will compare our results 
on the equation of state as a function of T /Tmin with the lattice data as a function 
of T/T,. 

The value of T^m is only dependent on the model parameter k in Eq.(^.l|). 
With given k = 0.43GeV, we can read the values of the critical temperature Tmin = 
201MeV and Tmin = 160MeV for c = + and c = — , respectively. Both values are in 



agreement with lattice QCD result on the critical temperature |2^, ^ 
4.2 The entropy density 




Figure 3: The scaled entropy density s/T^ as a function of scaled temperature T/Tmin 
with k = 0.43GeV and G^/L^ = 1.26 for c = + in (a) and c = — in (b), respectively. The 
dotted lines in (a) and (b) are lattice results as a function of T/Tc from p5[ |. 



Following the standard Bekenstein-Hawking formula [0, from the geometry 



given in Eq. ( p. 141) , one can easily read the black- hole entropy density s, which is 



defined by the area Aarea of the horizon: 

A 



'■area 




(4.19) 



Where is the Newton constant in 5D curved space and V3 is the volume of the 
spatial directions. It is noticed that the entropy density is closely related to the 
metric in the Einstein frame. 

With fixed k = 0.43GeV, Gg/^^ = 1-26, the scaled entropy density s/T^ as 
a function of scaled temperature T/Tmin is shown in Fig.^(a) for c = + (Tmin = 



-15- 



201MeV) and Fig.|(b) for c = - (T™„ = 160MeV), respectively. The dots in Fig.| 
are the lattice result taken from It is can be seen that the entropy density in 
both positive and negative quadratic correction hQCD models agrees well with the 
lattice result for pure SU (3) gauge theory. 



4.3 The pressure density, energy density and trace anomaly 

The pressure density p{T) can be calculated from the entropy density s{T) by solving 
the equation 

dp{T) 



dT 



sin 



(4.20) 





(a) 



(b) 



Figure 4: The scaled pressure density p/T^ as a function of scaled temperature T /Tmin 
with k = QA2>GeV and G^/L^ = 1.26 for c = + in (a) and c = — in (b), respectively. The 
dots are the lattice data as a function of T/Tc from |p5[. 



After integrating the Bekenstein-Hawking entropy density, the pressure density of 
the system can be obtained up to a integral constant Pq. In our numerical calculation, 
we have set po = to ensure piTmin) = 0. In some sense, this procedure is equivalent 
to choose a finite value of zir in Eq.( [4.17| ) to ensure AJ^ = at Tmin- 

The numerical result of the pressure density as a function of temperature is 
shown in Fig.|[ With fixed k = 0.43GeV, Gs/L^ = 1.26, the scaled pressure density 
p/T^ as a function of scaled temperature T/Tmin is shown in Fig.|(a) for c = + 
{Tmin = 201MeV) and Fig.|(b) for c = - [Tmin = 160MeV), respectively. The dots 
in Fig.||are the lattice result of pressure density for pure SU(3) gauge theory p5|. It 
can be seen that the pressure density in both hQCD models with c = ± agree well 
with the lattice result. 
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Figure 5: The scaled energy density e/T^ as a function of scaled temperature T/Tmin 
with k = OASGeV and G;^/L^ = 1.26 for c = + in (a) and c = —in (b), respectively. The 
dots are lattice data as a function of T/Tc from ||25|. 



Once we get the entropy density and pressure, we can get the energy density e 
immediately. The energy density is defined by 

e = -p + sT. (4.21) 

The numerical result of energy density as a function of the temperature is shown 
in Fig.^. The lattice data for pure SU{?>) gauge theory is shown with dots. It is 
found that the energy density from the two hQCD models with c = ± agrees with 
the lattice data. 

We also show the trace anomaly e — 3p in Fig.^. The trace anomaly shows a 
peak around T/T^m = 1-1 for the case of c = + and T/T^m = 1.2 for the case of 
c = — , and the height is around 3 for c = + and 2.5 for c = — . Both results are in 
agreement with lattice result. At very high temperature, the trace anomaly goes to 
zero, which indicates the system is asymptotically conformal at high temperature. 

4.4 The sound velocity and specific heat 

The sound velocity can be obtained from the temperature and entropy: 

c2 = ^^Qg^ = (4.22) 
^ d log s Tds I dT ' 

which can directly measure the conformality of the system. For conformal system, 
= 1/3, for non-conformal system, will deviate from 1/3. From Eq. (f4.22| ), we 



can see that the speed of the sound is independent of the normalization of the 5D 
Newton constant and the space volume V3. 
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• T 




(a) 



(b) 



Figure 6: The trace anomaly e — 3p as a function of scaled temperature T/Tmin with 
k = ^A2>GeV and G^jL^ = 1.26 for c = + in (a) and c = — in (b), respectively. The dots 
are the lattice data function of T/Tc from |25[]. 




(a) 



(b) 



Figure 7: The square of the sound velocity as a function of scaled temperature T/Tmin 
with k = 0.43Gey and G^/L^ = 1.26 for c = + {T^in = 201MeV) and c = - {Tmin = 
160MeV), respectively. The dotted line is the lattice data as a function of T/Tc from |p5|| . 



With known entropy density and the speed of sound, the specific heat C^, can be 
derived straightforwardly and takes the following expression: 



(4.23) 



At phase transition, the specific heat shows a clear A-type anomaly or divergent 
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(a) 



(b) 



Figure 8: The scaled specific heat C^jT'^ as a function of scaled temperature T/Tmin with 
k = 0.43Gey and G^/L^ = 1.26 for c = + (T™„ = 201MeV) and c = - {Tmin = 160MeV), 
respectively. 



behavior which is consistent with the jump or fast change of entropy density at 
Tc. Therefore, it can be used to determine the phase transition point. At very high 
temperature, the specific heat for the free gas of pure gauge SU{Nc) theory has the 
relation ||5^ 

a/T=^ = 4e/T^ ^ A{Nl - l)7rVl5, (4.24) 
and at large N^. limit [|], it takes the limit of 

CJT^ 4Ar27rVl5. (4.25) 

The numerical result of the square of the sound velocity is shown in Fig.^. At 
Tmm, the sound velocity square is around 0.05 for c = + and 0.02 for c = — , both 
are in agreement with lattice data 0.05. At high temperature, the sound velocity 
square goes to 1/3, which means that the system is asymptotically conformal. The 
numerical result of the specific heat is shown in Fig.|^. It can be clearly seen that 
the specific heat diverges at T^m. At T — i- oo, we can see that the scaled specific 
heat Cy/T^ approaches the free gas limit of pure gluon system 4(A''^ — l)7r^/15 with 
A^c = 3. 



4.5 Discussion 

In this section, we have investigated the equation of state for two hQCD models with 
positive/negative quadratic correction in the deformed warp factor of AdSs. It has 
been found that for both models, if we choose Tc = T^m, the properties of entropy 
density, pressure density, energy density and sound velocity are all in agreement with 
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lattice results for pure SU(3) gauge theory. As we have mentioned, Tc can be larger 
than Tmin- By comparing our numerical results with lattice data, it is observed that 
to fit the lattice data well, should be located in a neighborhood of Tmin from above, 
which means zjr should be in a neighborhood of zjji = 2.34GeV"^ from below. IfT, 
lies in this small region around Tmin, the numerical results on equation of state can fit 
lattice data almost equally well as that at T^ = T^in- It is worthy of mentioning that 
in our work choosing = T^in can fit lattice data better. Therefore, we identify 
Tmin as Tc in this work, and we believe that the qualitative results are not dependent 
on the exact location of T^. 

5. Heavy quark potential, Polyakov loop, spatial string ten- 
sion at finite temperature 

QCD vacuum is characterized by spontaneous chiral symmetry breaking and color 
confinement. The dynamical chiral symmetry breaking is due to a non- vanishing 
quark anti-quark condensate, (qq) ~ (250MeV)^ in the vacuum, which induces the 
presence of the light Nambu-Goldstone particles, the pions and kaons in the hadron 
spectrum. The confinement represents that only colorless states are observed in the 
spectrum, which is commonly described by the linearly rising potential between two 
heavy quarks at large distances, Ygqi^r) = ar, where r is the distance between quark 
anti-quark, and a (425MeV)^ is the string tension of the fiux tube. 

It is expected that chiral symmetry can be restored and color degrees of free- 
dom can be freed at high temperature and/or density. The chiral restoration and 
deconfinement phase transitions are characterized by the breaking and restoration of 
chiral and center symmetry, which are only well defined in two extreme quark mass 
limits, respectively. In the chiral limit when the current quark mass is zero m = 0, 
the chiral condensate (qq) is the order parameter for the chiral phase transition. 
When the current quark mass goes to infinity m — j- oo, QCD becomes pure gauge 
SU (3) theory, which is center symmetric in the vacuum, and the usually used order 
parameter is the Polyakov loop expectation value (L) [Q, which is related to the 
heavy quark free energy. 

In Sec.^, we have investigated thermodynamic properties of two deformed AdSs 
models with quadratic corrections, and have found that both models agree well with 
lattice results for pure SU (3) gauge theory. Therefore, we will not discuss the chiral 
phase transitions, but focus on deconfinement phase transition properties of these 
two models. In this section, we study the heavy quark potential, the Polyakov loop 
and the spatial Wilson loop at finite temperature, which are quantities related to 
deconfinement properties. It is worthy of mentioning that till now there is no good 
method to calculate the properties of these loop operators in the whole tempera- 
ture region within the framework of field theory, though chiral phase transition can 
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be described by using effective QCD models, e.g. the Nambu-Jona-Lasinio (NJL) 
model |55| and the linear sigma model [Q. The gauge/gravity duality offers a non- 



perturbative method to calculate the loop operators, in return, lattice QCD results 



[pYl , |58| , |59| of these loop operators will judge the validity of the hQCD model. 

At last, we want to mention that all the quantities of loop operator will be 
calculated in the string frame, which is different from the way of calculating the 
thermal quantities for equation of state. In order to avoid confusion and also to keep 
this section self-contained, we explicitly list the metric and the solved dual black-hole 
background from Sec.^, which will be used in this section: 







dsl 




As 


= cez\ 




pkz 




= l-fc 




Jo 



f{z)dt^ + + dx'dx' 



ixp ^2*^^ {Hc{x/k) — 1)J dx, (5.1) 

with c = ±, and /c and defined in Eq.(|3.3|) and (^.4]), respectively. 
5.1 The heavy quark potential 

We firstly study the heavy quark potential. The linear heavy quark anti-quark po- 
tential normally indicates the confinement of quarks in the vacuum. Above the 
deconfinement critical temperature T^, it is expected that the linear potential van- 
ishes and the Coulomb potential to be exponentially screened at large distance. The 
heavy quark potential at finite temperature has been analyzed in lattice QCD ||57 |. 



In the framework of gauge/gravity duality, we follow the standard procedure [pD 



and |6T| to derive the static heavy quark potential VqqIt) under the general metric 



background with a black hole, i.e, Eq. (|5.1|) . In SU{N) gauge theory, the interaction 
potential for infinity massive heavy quark antiquark is calculated from the Wilson 
loop 

W[C] = ^TrPexp[i j A^.dx'"], (5.2) 

where is the gauge field, the trace is over the fundamental representation, P 
stands for path ordering. C denotes a closed loop in space-time, which is a rectangle 
with one direction along the time direction of length T and the other space direction 
of length R. The Wilson loop describes the creation of a QQ pair with distance r at 
some time to = and the annihilation of this pair at time t = T. For T — )■ oo, the 
expectation value of the Wilson loop behaves as {W{C)) oc e~^^«'5. According to the 
holographic dictionary, the expectation value of the Wilson loop in four dimensions 
should be equal to the string partition function on the modified AdSs space, with 
the string world sheet ending on the contour C at the boundary of AdSs 

{W[C]) = Z!i,^^[C] ^ e-'--^^^ , (5.3) 
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where Sng is the classical world sheet Nambu-Goto action 

Sng = 77— / d^v^/DetXab, (5.4) 

with aq the string tension which has dimension of GeV^^, and Xab is the induced 
worldsheet metric with a, b the indices in the (77° = t,ri^ = x) coordinates on the 
worldsheet. 

We consider the following situation there are static quark and ant-quark linked 
by one string. The position of one quark is x = — | and the other is x = |. Under 



the background (|5.1|) , we can obtain the equation of motion: 



^ - Constant = v^—^, (5.5) 



. /l + M 



Here the r is dependent on zq which is the maximal value of z and z'{x = 0) = 0. In 
eq. (|5.5| ), we have defined fo = f{z = zq). For the configuration mentioned above and 
the given equation of motion, we impose the following boundary condtions z{x = 
0) = ^0,-2^(3; = ±f) = 0. Following the standard procedure, one can derive the 
interquark distance r as a function of zq 

r(zo) = [ du , ^^"'^ /.y^x -^TT-? — 7- (5-6) 

Jo 1 1 fo .A e^^^ f(zuiy)e^M-o'^) 

The heavy quark potential can be worked out from the Nambu-Goto string action: 



9q 



1 I ^2As{zou) 



^Zq Jo h _ _Jo_,Ae^±j^ V 

with Qq = It is noticed that the integral in Eq.( ^.7| ) in principle include some 
poles, which induces Vqq{z) 00. The infinite energy should be extracted through 
certain regularization procedure. The divergence of Vqq{z) is related to the vacuum 
energy for two static quarks. Generally speaking, the vacuum energy of two static 
quarks will be different in various background. In our latter calculations, we will use 
the regularized V^q' , which means the vacuum energy has been subtracted. 

Fig.|| shows the heavy quark potential Ygq/gq as a function of the distance 
between the two static quarks r for different temperatures T/Tmin = 1, 1.5,2,2.5,3. 
Fig.|^(a) is for c = +, i.e, the positive quadratic correction model, and Fig.|](b) is 
for c = — , i.e, the negative quadratic correction model. For c = +, the heavy 
quark potential has a linear part at zero temperature as shown in Ref. It is 

observed that when T > Tmin, the linear part vanishes, and with the increases of the 
temperature, the Coulomb part potential gets screened. This picture qualitatively 
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0.15 0.20 




(b) 



Figure 9: The heavy quark potential VQq/gq as a function of the separation distance r 
for different temperatures for c=+ (a) and c = — (b), respectively. Where k = 0.43GeV is 
used. 



agrees with lattice result in Ref. |57|. For the case of c = — , the heavy quark 



potential has no linear part at zero temperature as shown in Ref. |Q. It is seen 
that when T/Tmin > 1, the heavy quark potential only contains Coulomb part, with 
the increase of the temperature, the Coulomb part gets screened. 

We can see that, even though the heavy quark potential for hQCD models with 
positive and negative correction are quite different in the vacuum, the behavior of the 
heavy quark potential in these two hQCD models are qualitative the same above the 
critical temperature T^. The imaginary part of the heavy quark potential at finite 
temperature in the framework of AdS / CFT has been discussed in Ref. |p2l . 



5.2 The Polyakov loop 

The Polyakov loop is defined as 



L{T) = -tr Pexp 



l/T 



dt Ao 



(5. 



where N is the colors, the trace is evaluated in the fundamental representation and 
P stands for the path ordering. In QCD with infinitely heavy quarks, the Polyakov 
loop is related to the operator that generates a static quark |5^. We can interpret the 
logarithm of the expectation value ( L{T) ) as half of the free energy Fqq of a static 
quark-anti-quark pair at infinite distance. The Polyakov loop is an order parameter 
for center symmetry of the gauge group. In the confined phase, the 



QQ 



— )■ oo, 



which ensures {L(T) ) = 0, and thus the confined phase is center symmetric. The 
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deconfined phase is characterized by Fqq < oo and {L{T) ) ^ 0, which imphes the 
center symmetry is breaking in the ordered phase. 

In the framework of gauge/gravity duahty, the Polyakov loop is the Wilson 
loop wrapping the periodic imaginary time direction This operator should be 
computed in the string frame |6^. We follow the method in Ref. [0] to calculate 
the expectation value of the Polyakov loop, which is schematically given by the 
world-sheet path integral 

(L(r)) = / DXe-^^ , (5.9) 



where X denotes a set of world-sheet fields. is a world-sheet action. In principle, 
the integral ( ^.91) can be evaluated approximately in terms of minimal surfaces that 
obey the boundary conditions. The result is written as {L(T) ) = w„ exp[— S'„], 
where Sn means a renormalized minimal area whose weight is Wn- 

Given the background in the string frame Eq.( |5.1| ), we can calculate the expec- 
tation value of the Polyakov loop by using the Nambu-Goto action for 5*^ 



S 



1 



NG 



abi 



(5.10) 



with ap the string tension and Xah the induced worldsheet metric with a, h the indices 
in the {ij" = t,rf = z) coordinates on the worldsheet. It is noticed that here we have 



introduced a different string tension from that in Eq.( |5.4| ). This can be understood 
that along different direction, the string tension is different. We will also introduce 
another different string tension parameter for the spatial Wilson loop. This yields 



Q 9p 



Zh q2As 



dr ■ 



(5.11) 



where gp = A prime stands for a derivative with respect to z. 
The equation of motion for x takes the form of 



-f{z)x'/^i + fiz)ix'y 



0. 



(5.12) 



The above equation has a special solution x = const. Substituting this constant 



solution to Eq.( |5.1lD , one can get the minimal world-sheet. 



So = S'o + c. 



S' 



pi 



% / dz 



-,2As 



(5.13) 



Where Cp is the normalization constant which is dependent on the scheme of the reg 
ularization procedure. It is noticed that Eq.(|5.13|) takes the same form as that in 
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'g, = -o.i 



gp = -0.63 



(b) 



Figure 10: The expectation value of the Polyakov loop as a function of the temperature 
T/Tmin- (a) in the case of c = + with k = 0.43GeF, Cp = 0.1 and = -0.1,0.2,0.7, (b) 
in the case of c = - with k = QA'iGeV , Cp = 0.12, and gp = 0.2, -0.1, -0.63. The dots 



are lattice data as a function of T/Tc taken from |58|. 



where the black- hole background is taken from pure AdSs. Therefore, Performing 
the integral over z, we get the same analytic expression as in [|12 

gp ^/27[kzhErfi [\/2kzh) 



Sn 



vrT 
9p 



+, 



^kzhEri {V2kzh) 



ttT Zh 
Combining the weight factor with the normalization constant as Cp 
find 



(L(T))=exp{Cp-5;}. 



(5.14) 

(5.15) 
In Wo — Cp, we 

(5.16) 



The numerical result of the expectation value of the Polyakov loop as a function 
of the scaled temperature T/Tc is shown in Fig. [1^ (a) and (b) for positive and 
negative quadratic correction models, respectively. The dots are lattice data taken 
from Ref.|Q. For the positive quadratic correction model, i.e, c = +, the parameters 



k = 0.4:3GeV and Cp = 0.1 have been used for numerical calculation. We plot the 
Polyakov loop as a function of temperature for several values of gp = —0.1,0.2,0.7, 
it is found that gp = 0.7 can remarkably fit the lattice data for the pure SU{3) 
gauge theory. For the negative quadratic correction model, i.e. 



0.7 can remarkably fit the lattice data 

c = — , the pa- 
rameters k = OASGeV and Cp = 0.12 have been used for numerical calculation. 
We plot the Polyakov loop as a function of the temperature for several values of 
gp = 0.2, —0.1, —0.63, it is observed that gp = —0.63 can fit the lattice result quite 
well. 
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It can be seen that both models by choosing different parameters can agree well 
with the lattice result on the Polyakov loop. However, we have to point out that in 
the case of c = — , a negative value gp = —0.63 has been taken to fit the lattice data. 
This means the string tension ap introduced in Eq.( 5.1Cl|) is negative, which is not 



physical. Therefore, from the result of the expectation value of the Polyakov loop, 
the negative quadratic correction model can be excluded. Our result of the Polyakov 
loop for positive quadratic correction model fits lattice data remarkably. 

5.3 The spatial Wilson loop 

We now move on to discuss the spatial Wilson loop in the framework of the gauge/gravity 
duality. 

The non-Abelian gauge theories undergo a deconfining phase transition at high 
temperature. The physical string tension, characterizing the linear rise of the poten- 
tial between static quark sources with distance, decreases with increasing tempera- 
ture and vanishes above Tc. The potential becomes a Debye screened Coulomb po- 
tential in the high temperature phase, which is shown in Sec. p.l[ . However, the phase 
just above Tc is more complicated than a weakly coupled quark gluon gas because 
of the appearance of some non-perturbative soft modes, e.g, the Debye screening 
mass rrin ~ gT and the magnetic screening mass tum ~ g^T. The soft mode in 
magnetic sector cannot be handled in any perturbative scheme. It is found that the 
non-perturbative physics arising from the magnetic sector reflects the survival of the 
area law behavior for space-like Wilson loops, i.e. confinement of magnetic modes. 
The spatial string tension, ag, extracted from the area law behavior of spatial Wil- 



son loops, has been studied in lattice calculations and shows the expected 

dependence on the magnetic scale ^/o^ ~ g'^(T)T. 

In the framework of gauge/gravity duality, to calculate the spatial Wilson loop, 
we consider a rectangular loop C along two spatial directions (x, y) JlOl- The loop C 



should wrap around the which is a circle in the time direction. The expectation 
value of the loop can be calculated through the following AdS/CFT dictionary: 

K = y Me^^-. (5.17) 

Here X is the series of world sheet fields and stands for a world sheet action. 
Again, at large N, the saddle point approximation is valid, and we only need to 
take the minimum value of the Euclidean action among the saddle points. In fact, 
given the boundary conditions, there will be UV divergence from the field theory 
viewpoint. It is necessary to regularize and make it finite by a divergent subtraction 
by introducing a counterterm. 

Given the background in the string frame Eq. ( |5.1| ), we can calculate the expec- 
tation value of the spatial Wilson loop by using the Nambu-Goto action for Sw 



Sng = j d'^ri^/DetXab, (5.18) 
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with OL the string tension and Xah the induced worldsheet metric with a, h the indices 
in the [j]^ = x^tf = y) coordinates on the worldsheet. We take one of the spatial 
direction Y goes to infinity. The quark and anti-quark are set at Xj = ±r/2, we can 
get the Nambu-Goto action as 



r 

2 ^ h I z'"^ 



S='^Y dx-Jl + —- (5.19) 
2vr y_i z^\J f{z) 

here r stands for the separation between the static quark pair and = From 
the (|5.19|) , one can easily get the equation of motion for x: 

zz" + {f{z) + z'^) (2 - zd. In h) - ^z{z'fd. In f{z) = 0. (5.20) 



Following the standard procedure, one can obtain the simple equation: 

h ho 



(5.21) 



Where h{z) = e^^'^^^ and ho = h{zo). From the equation of motion, one can see that 

2z/2e"^^^''('''"^)^« 

r{zo) = zo =du, (5.22) 

•^0 ^Jf{zou) (1 - zy4g-4cfc2(,2_i),2^j 

where zq = z{x = 0) which defined by z'(a; = 0) = is the maximum value of z. The 
free energy can be read: 

SM = — d^- (5.23) 



After subtracting the UV divergence, one can get the area of the two dimen- 
sional minimal surface given by the classical configuration of the Nambu-Goto action 
Eq. (^.21|) , and the regularized spatial Wilson loop has the following form: 

Vs{zo) = ^[ t-A , ""^'^^^^^ =-lUz.-l|. (5.24) 



Fig. shows the regularized spatial potential Vg as a function of the distance 
between quark anti-quark r for several temperatures above T^m- For both c = + 
and c = — cases, the parameter k = 0.43GeV is used. It is found that for positive 
quadratic correction model, the spatial heavy quark potential always have a linear 
potential, and the slope of the linear potential increases with the temperature. This 
picture is in agreement with lattice result pR E^. For the case of c = — , the spatial 
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potential is flat at lower temperature, however, when the temperature is high enough, 
we observe that a linear potential rises up. 

Following the procedure introduced in Appendix y, we subtract the spatial string 
tension crs(T) and show the numerical results in Fig. |T^, where the lattice data for 
pure SU{2) and SU(?>) gauge theory are taken from |59| and [^| respectively. For 
the case of c = +, it is observed that the result of the spatial string tension fits 
very well with the lattice data in Ref.||25|[^. While for the model with negative 
quadratic correction, i.e. c = — , we find that the spatial string tension is not in 
agreement with the lattice date. 



6. Conclusion and discussion 

From the previous studies of constructing the holographic QCD models to describe 
hadron spectra and heavy quark potential, we have observed that a quadratic cor- 
rection is related to the confinement property of QCD. There are two methods to 
introduce the quadratic correction to hQCD models: the first method is to add it in 
the warp factor of the metric, the second method is in the dilaton background. How- 
ever, the connection between these two methods is lack of clarification. Moreover, 
the sign of the quadratic correction is still under debate. 

In this work, we have established a general framework for the minimal graviton- 
dilaton system. By solving the Einstein equations, we can self-consistently solve the 
dilaton background for any a given metric structure assumed in a phenomenological 
holographic model. Thus, the connection between the metric structure and the 
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Figure 12: The scaled spatial string tension ^fo's|Train as a function of the scaled temper- 
ature T I Train for the case of c = + in (a) with k = 0.43GeV and Qsv = 0.55 are used, and 
c = — in (b) with k = 0.43GeV and Qsv = 0.7 are used, respectively. The blue dots stands 
for lattice data for pure S'f7(3) gauge theory from [|^. The red dots are lattice data as a 



function of T /T^ for pure SU{2) gauge theory which is from [59| 



dilaton background is clearly revealed. Furthermore, we can observe: In the first 
method, the dilaton field should be solved out and its contribution must be taken 
into account; In the second method, the back-reaction of the dilaton background to 
the metric should be considered. In other words, both methods to introduce the 
quadratic correction are not self-consistent. 

In this work, we extend the graviton-dilaton system to finite temperature with 
a dual black-hole and self-consistently study the sign of the quadratic correction 
in the deformed warp factor by systematically investigating the thermodynamical 
properties and comparing with lattice QCD on the results of the equation of state, 
the heavy quark potential, the spatial Wilson loop and the Polyakov loop. We 
find that the bulk thermodynamical properties are not sensitive to the sign of the 
quadratic correction, and the results of both deformed holographic QCD models with 
a positive and negative quadratic correction agree well with lattice QCD results for 
pure SU(3) gauge theory. However, the results from loop operator favor a positive 
quadratic correction, which agree well with lattice results. Especially, the result from 
the Polyakov loop excludes the model with negative quadratic correction in the warp 
factor of AdSs. 

We would like to make one comment that the model with a quadratic correction 
in the warp factor is not equivalent to the model with a quadratic correction in 
the dilaton background. Therefore, our result does not indicate the failure of the 
KKSS model, where the negative quadratic correction is introduced in the dilaton 
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background. This deserves further study. 

It is interesting to ask the question why the bulk thermodynamical proper- 
ties, such as the energy density, sound velocity, are not sensitive to the sign of 
the quadratic correction in the deformed warp factor, while the thermodynamical 
properties of loop operators favor the positive quadratic correction model. The an- 
swer to this question is offered below. The crucial point lies in that fact that the 
bulk thermodynamical properties are defined in the Einstein frame, and the thermal 
properties of loop operators are defined in the string frame. In Fig.|l^ (b), we can 
observe that for positive and negative quadratic correction, the two warp factors in 
the Einstein frame are almost the same. Therefore, the bulk thermodynamical prop- 
erties are not sensitive to the sign of the quadratic correction. However, in Fig.|l3 
(a), we can observe that for positive and negative quadratic correction, the two warp 
factors in the string frame are quite sharp. This explains why loop operators are 
sensitive to the sign of the quadratic correction. 




c = - 




(a) 



(b) 



Figure 13: The warp factor in the string frame (a) and Einstein frame (b) for c = it. 
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APPENDIX 



A. Some simple analytical black hole solutions 

In this appendix, we work out two analytical gravity solutions of graviton-dilaton 
system by using Eq.( p.21| ) for given metric ansatz in Einstein frame Eq. (|2.14| ). Here 



we are interested in the solutions whose UV behavior is asymptotic AdS;^. We also 
impose the constrains which are /(O) = 1 near the z ~ 0, and requiring (f){z), f{z) to 
be regular at z = 0, Zh- 

To avoid misunderstanding, we should clarify that we will list two analytical 
gravity solutions in the Einstein frame. The metric ansatz in the Einstein takes the 
following form: 

2 ^2g2A,-M / 2 dz"^ 



dsj, = \^-fiz)dt' + + dx'dx'j (A.l) 

-f{z)de + ^ + dx'dx' . (A.2) 



V ^ ' f{z) 

Where As{z) and 0(z) is defined by Eq.(|21^) and Eq. (|27[3| ) and Ae{z) = As{z) - 

3 ■ 

Using the general solutions in Eq.( p.2l| ), one can easily find some exact new 
gravity solutions. Here, we do not repeat the procedure and details of the calculation 
to obtain the new gravity solutions, we just list two simple analytical gravity solutions 
for example. 

The first solution is given as following: 



Ae{z) =\oi 

yzosmh(-j J 
f{z) = 1-^(2-3 cosh(-) + cosh3(-)), 

O Zq Zq 

m = 1^- (A.3) 

AZq 

Where zq is integral constants and Vi is the constant from the dilaton potential 

Wi(sinh^(|)) + 9(sinh^(f)) + 12 
Ve(0) = ^ j-^ ^ . (A.4) 



The solution above is given by the graviton-dilaton system (|2.12 ) with nontrivial 



dilaton potential Ve- Where the subscript E stands for the dilaton potential in the 
Einstein frame. Here we should emphasize that the exact relation of the dilaton 
potential between the Einstein frame and string frame is given by the Eq. ( |2.11| ). 
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The second new gravity solution is as below: 

Ae{z) = -log + (A.5) 

/ 7 6 o 5 4 \ 

= 3y2sinh~^ (V?) ■ ^'^■'^^ 
Where zq is integral constants and V2 are the constants from the dilaton potential 



Ve{<P) 



3V-.sinh^^(^) 3V-,sinh^°(^) 
35L2 10L2 10L2 

42sinh^(^) 42sinh^(5^) ^2 



L2 L2 L2 ■ 

In this context, we do not discuss the application of the above two gravity solu- 
tions, such as the thermal dynamical quantities and stability of the two black hole 
solutions, which are left to the future work. 



B. Free energy density at UV boundary for the hQCD model 
with As{z) = c/c^z^ 

We offer the details to derive the free-energy density at UV boundary z = e ^ for 
the hQCD model with ^^(,2) = ck'^z'^ described in Section ^. 

For black-hole case, the solutions in Section ^ expanded around UV boundary 
have the following expressions: 

0(.) = + ^k'z' + ^k'z' + ^kh' + o{z% (B.l) 

f{z)= l-ilk^z^-^k'z' + o{z^), (B.2) 

AJz) = --k'z^ - —k^z"^ - — A;«^« + o(z^). (B.3) 
^ ^ 5 105 315 ^ ' ^ ' 

For the case of thermal gas, the solutions at UV boundary take the expression 

as: 

.2 2 ^ p|g^8-105i^) 6 6 

(Poiz) = P2k z H k z H k z 

^ ^ 15 3780 

Uz)= 1, (B.5) 

^ 45 2835 25515 ^ v ; y J 



- 32 - 



To compare the free-energy density of the black-hole case and the thermal gas 



case, one has to impose the following conditions |^ : 



0o(^) |e= m U (B.7) 
Pho{z) 1= Ph{z) |„ (B.8) 
V^hl{z) 1= V,h\z) I . (B.9) 

Comparing Eg. ( p.3| ) with Eq. (p.6| ), we can get the relationship between p2 and k, 

F2 2 y IQ 8960 ) 

Then we get the free energy densities for the black-hole and thermal gas at UV 
boundary as follows: 

= 2M^ {?,h{eff{e)h\e) + \h{eff{e)^ 

= _ Qa;U3^3 _ Ik'L^M^fl^ + o(e2) (B.ll) 



TG 



6L^M^ /6,.^,,,, 3 



(^h^L^M^ - h'L^M^f^^ + o{e'). (B.12) 



From the above result, one can derive the free energy density difference as follows: 
= -h^M'L'f^ - 2M'blizjnMzrn). (B.13) 



From Eqs.(p.3|),( [4.3| ) and ( [4. ID]) , one can express the temperature dependent coeffi- 
cient with the entropy density and temperature as 

As a result, one gets 

= ~Ts - 2M%{zjn)b',{zjR). (B.15) 
C. To extract the spatial string tension 

In this appendix, we show how to extract the spatial string tension. A fist-principle 
method for determining the spatial string tension in hot QCD matter is to measure 
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the rectangular spatial Wilson loops as we show in Sec. 5^. The spatial potential 
Vs{r) is defined by: 

Vs(r) = - lim -\ogW(r,Y). (C.l) 

y-s-oo r 

The spatial string tension can be obtained: 

a,(T) = lim (C.2) 

r— !>oo r 

To extract the spatial string tension, one can use the following relation: 

Vs = S,y r^asT r > 1. (C.3) 
Under the general background 

ds^ = -At- f{z)dt^ + TTT + dx'dx' , (C.4) 
^ V f{z) J 

the separation of the two static quarks has the form of: 

rizn) = zn I — , du. (C.5) 

Jo y/f {uZo) (1 - zy4g4A,(.o)-4A.(..o)) 

The free energy of the spatial Wilson loop is given by: 

gsv 1 ^2As(uzo) 

Sxvizr,) = / — — du. (C.6) 

nzo Jo Z/2 (lyzo) (1 - i.4e4A.(.o)-4A.(..o)) ^ ' 

We assume that the function As{zq) is a regular function here. One should focus on 
the denominator of the Eq. (|C.5| ) and Eq. ( p.6|) which make main contribution to the 
integral when z/ ~ 1. If one just focus on the dominant part, one can expand the 
denominators at z/ = 1 respectively. After the expansion up to the order 0((z/ — 1)^) 
and integral out the formula, one can find the coefficient: 

(l-z/)(-4^o/(^o)^;(^o) + 4/ (^o)) 

+ (1 - z/)' (4^oV' (^o) K (^o) - 8^oV (^o) A'^ (^o) ' + 24^0/ (^^o) A'^ (^o) 
+2z'j (zo) A': (zo) - Azof (zo) - Uf (zo) ) (C.7) 
In term of the following useful integral formula: 



1 _ 2ArcTanh 



v^a(l - z/) + 6(1 - z/)2 Vb 



a > and 6 > 0, (C. 



one can notice that the integral will be divergent when a = 0. r{zQ) can be obtained 
and has the form of 



Azq tanh ^ ^- 



r(.o) = '-^^ + 0(1) (C.9) 
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with cl and 6^ the coefficients appeared in Eg. ( |C.8|) respectively. 
Following the same strategy, we can get the spatial potential 

7TZQy/0v 



Where ay and by are the coefficients mentioned in Eq.( |(J.(j| ) respectively. Comparing 
Eq. (|C.9D and Eq. (|C.10|) , one can easily get the spatial string tension. Up to the order 



0(1 — z/)^, the coefficients listed in Eq. |C.7| have the following exact relations such as 
ctL = ay and = by in Eq. (|C.9|) and Eq. (|C.10|) . The spatial string tension has the 



form of 

cTs/gsv = — — (C.ii) 

Here, One can find that the spatial string tension as depends on the warp factor 
As{z). In pure AdS^ case, 

(^s/9sv = (C.12) 
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